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Abstract 

Faddeev equations for elastic Nd scattering have been solved using modern 
NN forces combined with the Tucson-Melbourne two-pion exchange three- 
nucleon force, with a modification thereof closer to chiral symmetry and the 
Urbana IX three-nucleon force. Theoretical predictions for the differential 
cross section and several spin observables using NN forces only and NN forces 
combined with three-nucleon force models are compared to each other and to 
the existing data. A wide range of energies from 3 to 200 MeV is covered. 
Especially at the higher energies striking three-nucleon force effects are found, 
some of which are supported by the still rare set of data, some are in conflict 
with data and thus very likely point to defects in those three-nucleon force 
models. 
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I. INTRODUCTION 



One major goal in nuclear physics is to establish the properties of nuclear forces and to 
understand nuclear phenomena by solving the many- nucleon Schrodinger equation driven by 
those elementary nuclear forces. Meson theory had an important impact for the construction 
of nuclear forces (both NN and 3N forces), but it lacks systematics like it would be given by 
an expansion parameter and the meson-nucleon vertices had to be parametrised in an ad-hoc 
manner. Nevertheless the one-pion exchange is undisputed and the most advanced formu- 
lation of meson-exchanges in the so called full Bonn potential [|TJ is remarkably successful. 
Because of its energy dependence - a consequence of deriving it by old fashioned time or- 
dered perturbation theory - it is not useful in A > 2 systems. Energy independent one-boson 
exchange versions thereof, however, are useful |2|U . In addition more phenomenological NN 
potentials have been constructed with the aim to describe the rich set of experimental NN 
data as precisely as possible. This leads to an often called new generation of realistic NN 
potentials: AV18 0, CD Bonn §, Nijm I, II and 93 §. They describe the NN data set 
with an unprecedented precision of x 2 P er data point very close to one. Very recently an 
updated CD Bonn |J appeared, which takes newest data into account but has not been used 
in the present article. An upcoming approach to construct nuclear forces in a systematic 
manner is chiral perturbation theory |7]~ni|. First applications to three- and four-nucleon 
systems have been done |14] . 



In recent years it became possible to solve exactly three- and four-nucleon bound states 



using standard integration and differentiation methods [15,16|. Stochastic techniques allow 



to go beyond A=4 and right now low energy states of nuclei up to A=8 are under con- 
trol [T7|,18|. In all cases those realistic NN forces fail to provide the experimental binding 



energies; there is clear cut underbinding. For instance this amounts to 0.5-1 MeV in case of 
three nucleons and to 2-4 MeV in case of 4 He. A natural further step is the consideration 
of 3N forces. This is an even harder theoretical challenge and presently the most often used 
dynamical process is the n — n exchange between three nucleons with an intermediate excited 



nucleon state, the A ||19|| . This is augmented by further ingredients of various types as will 
be detailed below. By properly adjusting parameters one can achieve correct 3N and 4N 
binding energies and reaches even a fairly successful description of low energy bound states 
energies of up to A= 8 [2~0"|. However, in the latter case the results point to an insufficient 
spin-orbit splitting of nuclear levels in light nuclei as e.g. in 5 He JT7|. This may be caused 
by a wrong spin structure of present day 3NF's or by not well enough established 3 Pj NN 
force components. It will be interesting to see in the future the predictions of nuclear forces 
based on chiral perturbation theory. 

Though this first signal on 3N force effects resulting from discrete states is important, a 
more detailed investigation of 3NF properties can be carried through in scattering processes, 
where a rich set of spin observables is available. The tremendous advance in computational 
resources allowed in recent years to make exact predictions for three- nucleon scattering using 
nuclear forces in all their complexities [21]. Also experimentally one can access nowadays 
spin observables in Nd scattering where in the initial states the deuteron and / or the 
nucleon is polarized and after the reaction also the polarization of the outgoing particles 
can be measured p2] , |24] -pQ| . This leads to a very rich spectrum of observables in Nd elastic 



scattering and the Nd break-up processes. Such a set of spin observables will be a solid basis 
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to test the 3N Hamiltonian. Using available model Hamiltonians one can provide guidance in 
selecting specific observables and energies which are most appropriate to see 3NF properties. 
It is the aim of this article to do exactly that and to compare the theoretical predictions 
with already existing data. 

In section [Tl] we review briefly our 3N scattering formalism and display the 3NF model 
forces which are presently en vogue and which we use. Some technical details referring to 
the partial wave decomposition of the momentum space representation of the Urbana IX 
3NF |31j are given in the Appendix. In this article we restrict ourselves to elastic Nd 



scattering and refer to a forthcoming article for the break-up process. Our predictions 
for various nuclear force combinations and the comparison to available data are given in 
section ffll. We conclude in section ffV. 



II. SCATTERING FORMALISM AND 3NF MODELS 



21 



for a general overview on 3N scattering and specifically our way to 

It 



We refer to 

formulate it. For the inclusion of 3NFs we found meanwhile a more efficient way [|32 



is a direct generalization of what is being used for the 3N bound state ||33 
amplitude T via our central equation 



We define an 



T = tP4> + (1 + tG ) V} 1] (1 + P) <p + tPG Q T + (l + tG )V^> (1 + P)G T (1) 

The initial channel state occurring in the driving terms is composed of a deuteron and 
a momentum eigenstate of the projectile nucleon. The NN t-operator is denoted by t, the 
free 3N propagator by Go and P is the sum of a cyclical and an anticyclical permutation of 
three particles. The 3N force V4 can always be decomposed into a sum of three parts 



(i) 



v 4 = vP + V, 



r(2) 



v.. 



(3) 



(2) 



where V4 is symmetrical under the exchange of the nucleons jk with j ^ i ^ k. As seen 
in Eq. (|l|) only one of the three parts occurs explicitely, the others via the permutations 
contained in P. The physical break-up amplitude is given via 



U = (1 + P)T. 



(3) 



The Faddeev-like integral equation (HD has the nice property that its iteration inserted 
into Eq. (|3|) yields immediately the multiple scattering series, which gives a transparent 
insight into the reaction mechanism. Here in this article we concentrate on elastic scattering, 
whose amplitude is given by 



U 



PGo 1 



PT + V 4 (1) (1 + P) <p + V^> (1 + P) G T 



(i) 



(4) 



The first term is the well known single particle exchange diagram, then there are terms where 
either V4 or the t's interact once and then the remaining parts result from rescattering among 
the three particles. Again inserting the iteration of T as given in Eq. (|IJ) into Eq. (^) yields 
a transparent insight [3~I . 

The definition of the various spin observables can be found in [21,35|. They have the 
general form 
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( 5 f 1 = \ E( s ")^ (MFMiS") , (5) 

where I is the elastic cross section summed over the spin states in the final state, M is the 
physical elastic scattering amplitude related directly to U and is a suitable set of 3N spin 
operators. 

We shall encounter nucleon and deuteron vector analyzing powers A y (N) and A y (d) 
(iTn), where in the initial state either the nucleon (N) or the deuteron (d) is vector polarized. 
Further, the deuteron can be tensor polarized in the initial state leading to the three tensor 
analyzing powers T 2 fc (k = 0,1,2). Also both particles can be polarized in the initial state 
leading to very many spin correlation coefficients C a ^, where a refers to the spin directions of 
the nucleon and beta to vector and tensor polarizations of the deuteron. Further information 
on the dynamics can be found in spin transfer coefficients K@ , where a describes either a 
polarized nucleon or a polarized deuteron in the initial state and /3' similarly the polarization 
for a particle in the final state. Of course all those quantities depend on the scattering angle. 

Our nuclear model interaction consists of one of the NN forces mentioned in the intro- 
duction and a 3NF. For the 3NF we use the 27r-exchange Tucson Melbourne (TM) model, 



a modified version thereof and the Urbana IX force. The TM model |3b] has been around 
for quite some time. It is based on a low momentum expansion of the tt — N off (the- mass) 
shell scattering amplitude. It has the form 



r 2 • f 3 ( a + bQ ■ Q' + c(Q 2 + <?')) + di f 3 x f 2 • n a x ■ Q x Q' 



(6) 



The elements of the underlying Feynman diagram are obvious: the two pion propagators 
depending on the pion momenta Q and Q', the two tcNN vertex amplitudes and most 
importantly the parametrisation of the irN amplitude inside the curly bracket which is 
combined with the isospins r 2 and r 3 of the two accompanying nucleons. On top of all that 
there is a strong form factor parametrisation given by 



h(q- 



'A 2 



mz 



A 2 + Q : 



(7) 



In what we denote by the TM 3NF we use the original parameters a = 1.13/m,,-, b = 
— 2.58/m^, c = 1.0/m^, d = — 0.753/m^. They incorporate among others the physics 
resulting from an intermediate A in a static approximation. The cut-off parameter A is used 
to adjust the 3 H binding energy separately for different NN forces |37| . For the convenience 
of the reader we show the A- values in Table 1. Of course in a meson exchange picture 
additional processes should be added containing other meson exchanges like tt — p, p — p; 
also different intermediate excited states might play a role. To some extent 3NF models 
with respect to those extensions have already been developed and applied |38|-44l|. Further 
studies should be performed. 

The parametrisation of the TM 3NF has been criticized somewhat, since it violates 



chiral symmetry |42|,|43[]. A form consistent with chiral symmetry (though not a complete 
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one to that order in the appropriate power counting) is obtained by modifying the c-term 
so that the long-range part is absorbed into the a-term, leading to a new a' = a — 2m\c = 
— 0.87/m^ [pE2| , pE3"| what essentially means a change of sign for a and that the short range 
part is dropped. This form will be called TM' later on. The corresponding A-value when 
used with the CD Bonn potential is A = 4.593 m?,-. 



The two-meson exchange 3NF has also been studied by Robilotta et al. ||4 leading to 
the Brazilian 3NF. It is similar to the one of TM and also the results gained for low energy 



Nd elastic scattering observables are similar to the ones for the TM 3NF. In this article 
we do not take that force into account. Instead we included the Urbana IX 3NF pi] , which 
is heavily used in the Urbana-Argonne collaboration. It will be interesting to see its effects 
for 3N scattering observables. At very low energies it has been used in that context before 
by the Pisa group [£5|. That force is based on the old Fujita-Mijazawa ansatz [|7|| of an 



intermediate A occurring in the two-pion exchange and augmented by a spin independent 
short range piece. It has the form 



Va 



(i) 



A 



2 IT 



{X 12 , X 13 } {fi • r 2 , T\ ■ r 3 } + - [X12, X 13 ] [t x ■ f 2 , n ■ f 3 ] 



UoT*(r 12 )T*(r 13 ), 



(8) 



where 



(9) 



with 



YJr) 



(l-e- cr 



(10) 



and 



TJr) 



1 + 



■m w r (m^r) 2 



(11) 



and where 



3 a. 



— <Ji 



(J,; 



(12) 



is the tensor operator. 

Since we work in momentum space and in a partial wave expansion the form given in 
Eq. (H) has to be rewritten. We could follow the steps laid our before for the corresponding 
representation of the TM 3NF [48[] and delegate all that to the Appendix. 

Since there is no apparent consistency of the mostly phenomenological realistic NN forces 
and the 3NF models we test various combinations thereof. In all cases, however, we require 
that the particular choice for the 2N interaction and the 3NF should reproduce the exper- 
imental triton binding energy. Some of the 3N observables scale with the triton binding 



energy |49|]. The adjustment to the triton binding energy has the advantage that our inves- 
tigation is not misled by these scaling effects. 
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With respect to the intermediate A one should say that very likely the static approxi- 
mation is not justified and the A should be allowed to propagate like the nucleon. This has 
been pursued intensively for instance by the Hannover group |5(| and recent work has been 
also devoted to the 3N continuum []5TJ . 

In view of all that it is quite clear that our present study is not at all complete but can 
at most provide some insight, what kind of effects specific 3NF models might generate. As 
we shall see effects of that sort are needed, since NN force only predictions often fail to 
describe the data, especially at the higher energies. These challenges call for a systematic 
approach and at least for the leading spin structures chiral perturbation theory might be a 
good candidate J8| . This is left to a future investigation. Here we concentrate on the current 
models and show their strengths and failures. 



III. PREDICTIONS OF 3NF EFFECTS AND COMPARISON TO DATA 

Since we would like to cover a wide range of incoming neutron energies from below the 
nd breakup threshold up to 200 MeV it is necessary to take a sufficient number of partial 
wave states into account in order to get converged solutions of the Faddeev equations. In all 
calculations presented in this paper we went up to the two-nucleon subsystem total angular 
momentum j max = 5. This corresponds to a maximal number of 142 partial wave states 
(often called channels) in the 3N system. We checked that the convergence has been achieved 
by looking at the results obtained for j max = 6, which increases the number of channels to 
194. This convergence check refers to a calculation without a 3NF. The inclusion of 3NF's 
has been carried through for all total angular momenta of the 3N system up to J = 13/2. 
These high angular momenta are required at the higher energies > 100 MeV. The longer 
ranged 2N interactions require states up to J = 25/2 at the higher energies in order to get 
converged results. 

A phenomenological criterium for 3NF effects is that the data lie outside the spread of 
NN force predictions only. In the following figures we shall always include a shaded band 
(called "bandl"), which covers the predictions of the AV18, CD-Bonn, Nijml, II and 93 
NN forces. Unfortunately we cannot include the pp-Coulomb force in our approach and 
thus have to live with some theoretical uncertainty when comparing to pd-data. At the 
higher energies, however, those effects should be small. Also in case of AV18 we do not take 
the various electromagnetic corrections into account, which leads for example to a slightly 
wrong deuteron binding energy (E^ = 2.242 instead of 2.225 MeV). This, however, has only 
a small effect on our results, which is mostly of kinematical origin, since the phase-shifts 
obtained without those additional terms differ only slightly from the standard ones. The 
kinematical effects are seen predominantly in the breakup process, where this small defect 
in the deuteron binding energy leads to correspondingly small shifts in peak structures like 
for instance final state interaction peaks. 

We shall combine another group of curves into a band. The TM 3NF can be combined 
with the five NN forces. In all cases the cut-off value A in Eq. (0) has been adjusted 



separately for each NN force to the H binding energy ||37|| . Since that interplay is a purely 
phenomenological step the outcome is theoretically not under control and we combine all the 
results into a second band (called "band2" in the following). Next we want to compare the 
TM 3NF and the modified TM', which is more consistent with chiral symmetry. We combine 
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it with CD-Bonn and show the CD-Bonn+TM' prediction as a dashed curve. Finally we 
compare the TM and the Urbana IX 3NF's and combine them with AV18. The combination 
AV18+URBANA IX will appear as a solid curve. There are clearly more combinations 
possible but it is sufficient to get an orientation on the magnitudes of expected effects. 

We begin with the differential cross section in Fig. [I]. The various NN force predictions 
are rather close together with a small spread in the minima. Including the TM 3NF there 
is again a small spread in the minima (practically negligible at 3 MeV) but the minima are 
shifted upwards, rather well into the data except at 3 MeV, where the 3NF prediction 
is shifted slightly downwards. The phenomenon of shrinkage of the spread between different 
NN potential predictions by including a 3NF is often called a scaling phenomenon. It occurs 
at low energies and is related to the three-nucleon binding energy, which by construction 
is common to all those curves in band2. The TM and TM' 3NF's together with CD-Bonn 
give slightly different predictions in the minima especially at the two highest energies. (The 
CD Bonn+TM prediction lying inside band2 is not shown.) In the backward angular region 
they differ significantly and the 135 MeV precise backward angular distributions data prefer 
the TM 3NF. (See insertion in Fig. 1; again the CD Bonn+TM prediction is not explicitely 
shown.) On the other hand TM and URBANA IX together with AV18 are very similar 
at the two higher energies but differ significantly at E— 65 MeV. Certainly Coulomb force 
effects should be taken into account at this energy before a final conclusion can be drawn. 
The few nd data near the minimum would strongly disagree with all our 3NF predictions 
and a confirmation (or rejection) would be highly desirable. However, independent from 
possible Coulomb force contributions, it is clearly seen, that even such a simple observable 
as the elastic scattering differential cross section exhibits large 3NF effects at higher energies. 
These effects are not trivial and depend not only on the incoming energy and the angle but 
also on the particular 3NF used. This calls for precise data for this observable in order to 
study the 3NF properties. 

Let us now regard a selection out of the many spin observables in elastic Nd scattering. 
Fig.|2] shows A y (N). The band for NN force predictions is always rather narrow, whereas 
band2 for the lowest and highest energy is distinctly broader. The two bands are separated 
predicting clearly 3NF effects especially at higher energies. The TM and TM' predictions 
are distinctly different as well as the TM and Urbana IX predictions. It is interesting to 
note that here TM' with CD Bonn and Urbana IX with AV18 are very similar (except at 
the lowest energy) and predict only small 3NF effects at 65 MeV, which are compatible with 
the A y data at this energy. At higher energies their effects become quite different from the 
TM ones. While in the region of the A y minimum around 8 cm ~ 110° they increase A y as 
compared to the pure 2N force predictions, their action decreases A y in the backward angular 
region, contrary to the action of the TM 3NF model, bringing the theory closer to the data. 
In this way the TM' and Urbana IX 3NF's seem to solve partially the A y problem found 



at higher energies in P7LE3I. At 3 MeV the clear discrepancy of all theoretical predictions 



to the nd data of [54 is seen. At such a low energy it is well known that Coulomb force 



effects are large for A y decreasing significantly its maximum when compared to nd data [5"o* 
Thus also pd data lie very clearly (due to much smaller error bars) above all theoretical 
predictions. We see that this very well known low energy A y puzzle [|5jj cannot be solved by 
the 3NF models we are using in this article. A slightly increased maximum of A y for TM' is 
far too small to play any significant role and possibly the solution should be also sought in an 
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improvement of the 3 Pj NN force components |57] to which low energy A y is very sensitive. 
We would also like to point to the very recent result based on chiral perturbation theory 
HI , where A y can be described quite well in next-to-leading order (NLO). In that order of 



the power counting 3NF's do not yet contribute. Those effective chiral NN forces are very 
different from the conventional ones. But this NLO result is just an intermediate step and 
the final answer has to wait for higher order contributions, which improve systematically 
the observables in the 2N, 3N, ... systems at the same time. 

For iT\\ shown in Fig. || the two bands are distinctly different and clearly the TM-band 
is supported by the data at the higher energies. In that case TM and TM' are close together 
and also TM and URBANA IX except around 120 at the highest energy. The data shown 
in the figure for 190 MeV are taken at 197 MeV. Unfortunately they are absent around 
120 °. 

Next we regard the three tensor analyzing powers T 20 , T 21 and T 22 in Figs. J§-|| For T 20 
the situation is very challenging. At 135 MeV the data between about 120 °-150 ° do not 
agree with the overlapping bands and above 150 ° they lie just between the two bands. At 
small angles they agree with the overlapping bands and follow then the NN force prediction. 
In addition at the higher energies TM and TM' differ as well as TM and URBANA IX. The 
strong deviations of the theory to the data at 3 MeV is simply caused by Coulomb force 
effects H . 



For T 21 the situation is different. The two bands are clearly distinct. Again the different 
3NF's predictions deviate strongly from each other. The data at 135 MeV follow more bandl 
than band2 and at the small angles TM' or URBANA IX are preferred. Clearly this is a 
rather contradictory situation. 

For T 22 the bands differ but the special 3NF predictions (dashed and solid lines) are 
similar but lie outside band2. Except at very backward and forward angles, where all curves 
essentially coincide, there is disagreement with the data at E— 135 MeV. For all three T 2 fc's 
data at 190 MeV would be very valuable, since the various theoretical predictions differ 
dramatically. 

There are many spin-transfer coefficients and we selected more or less arbitrarily five of 
them. In Figs. |T|-|TT| we show Ky Z , K^ y) K%' x , K^ z and K%' z ' . For Ky Z the bands strongly 
deviate at the higher energies and TM and TM' as well as TM and URBANA IX differ 
drastically. The deviation of the bands from each other is less pronounced for and also 
the different 3NF predictions are less distinct except at the highest energy. Two data points 
agree with the 3NF predictions, while one, at 150 °, is below any theoretical prediction. For 
K% x the bands differ only at the lowest energy (E= 3 MeV) significantly. Otherwise the 
3NF effects are rather modest and also differences among the different 3NF's. The data at 
135 MeV go across the various predictions. For K% z the bands differ at the high energies 
and the two special 3NF predictions deviate significantly from the TM predictions. The one 
data point at 135 MeV lies somewhere in between. Finally K% z show again dramatic effects 
in relation to the two different bands. Also the special 3NF predictions are clearly different 
at the two higher energies. 

Lastly we regard in Figs. |l2| - p^ three different spin correlation coefficients C xy>x , C yy and 
C zz . The effects are dramatic for C xy ^ x : the bands are quite different and also the special 
3NF predictions. For C yy there are data |2(| at E p = 197 MeV, which we inserted into the 
figure for 190 MeV. The data support the curves inside band2. The AV18+URBANA IX 
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is significantly closer to experiment than CD Bonn+TM'. Finally for C zz the bands differ 
a lot at the high energies and the special 3NF predictions differ from the TM ones as well. 
One notes that URBANA IX and TM' are similar. 

IV. SUMMARY AND CONCLUSIONS 

We performed a study of popular present day 3NF models with respect to the effects they 
cause in the 3N continuum. Based on the comparison of the realistic NN force predictions 
alone ("bandl") to the predictions of all NN forces combined with the TM 3NF ("band2") 
one sees in many spin observables very drastic effects, which should clearly be discernible 
by experiments. On top the three different 3NF models, TM, TM' and URBANA IX com- 
bined (arbitrarily) with CD-Bonn and AV18 lead again to other very distinct predictions. 
Specifically the effects are angular and energy dependent becoming large at higher energies. 
It seems that with a sufficiently rich and precise data basis such diversity of effects should 
allow to nail down the proper spin structure of 3NF's. 

Unfortunately there are up to now only few data available. The ones for the differential 
cross section support the shift in theory caused by 3NF's. The existing high energy cross 
section data in the backward angular region prefer the structure of the TM 3NF. Also 
the existing deuteron vector analyzing powers at higher energies support rather well the 
predicted 3NF effects. On the other hand this three-body interaction predicts too large 
effects for the nucleon analyzing power. This observable seems to prefer the modified version 
of the Tucson-Melbourne model, TM', which is consistent with chiral symmetry or the 
URBANA IX. For the tensor analyzing powers the situation is totally chaotic, for some 
scenarios one finds agreement, for others a strong disagreement: there is no preference for 
any of them. Clearly we are at the very beginning in investigating the spin-structure of the 
3NF. Nevertheless the effects of all those 3NF's are typically of the right order in magnitude, 
when they can be checked against data but the signs are not yet under control. The spin 
transfer coefficients carry also a lot of information and in some of them the two bands differ 
very much. Finally, the spin correlation coefficients appear also very informative and the 
very first data for C yy support the TM and URBANA IX. 

Altogether the only conclusion possible is that the most popular current 3NF models 
show a lot of effects and data are needed to provide constraints. There is hope that further 
theoretical work guided by the chiral effective field theory approach will help to establish 
the proper spin structure of the three-nucleon force. 
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APPENDIX A: PARTIAL WAVE DECOMPOSITION OF THE URBANA IX 3NF 

IN MOMENTUM SPACE. 



The Urbana 3NF in Eq. @ has to be put into a form suitable for the evaluation in partial 
wave decomposition. Therefore we rewrite Eq.(|8|) using the (anti) commutator of the isospin 
operators 

{fi ■ f 2 ,fi • t 3 } = 2 r 2 • r 3 
[fi • r 2 , T\ ■ f 3 ] = 2i n ■ f 2 x r 3 (Al) 

V4 ^ reads then 

V 4 (1) = A 2n 2^t 2 -t 3 + ^t 1 -t 2 x r 3 ) X 12 X 13 

+A 27r 2^-7^-^fi-7SX7|J X i3 Xi 2 

+C/ T 2 (r 12 ) T 2 (r 13 ) 

= A 2t 2 ^ • 73 + ^ 71 • f 2 X 7^ Xi 2 Xi 3 
+A 2?r 2^-^-^fI-f^X7^j Xi 3 Xi 2 

+y (^ 2 (r 12 ) T w 2 (r 13 ) + T 2 (r 13 ) T 2 (r 12 )) (A2) 

Because the T 2 operators commute which each other, we could choose the symmetrical form 
in the last line. 



The aim is to find matrix elements with respect to our standard basis states [59 



\p qa >i= \p q {ls)j (A^) J (jJ)JM {t~)TM T > t (A3) 

Here the magnitudes of the Jacobi momenta of the subsystem and the outer particle are p 
and q, respectively. The angular dependence is expanded in partial waves. Corresponding 
to p and q we introduce angular momenta / and A. As indicated in Eq . (pOl) the orbital 
angular momenta couple with the spin of the subsystem s and the outer particle \ to the 
total spin of the subsystem j and outer particle J. These angular momenta are combined to 
the total spin J and its third component M.. The total isospin of the subsystem t couples 
with the isospin of the outer particle | to T and its third component M T . Because we will 
make use of several sets of Jacobi momenta, we append an index i which gives the number 
of the outer particle. 

According to Eqs.(|l|) and (|4|) the operator V4 acts on completely antisymmetric 3N 
states |x >= (1 + P) \<f) > or \x >= (1 + P) G T. In the next paragraphs we would like to 
establish some consequences of this fact. 

To that aim we introduce a short hand notation for our basis states: 

\(Jk)i>= \p qa>i (A4) 

In contrast to Eq.( [A3| ) this definition fixes the ordering within the subsystem, (ij) is un- 
derstood to fix the momentum vector to p = \{ki — kj) and the spin and isospin coupling 
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within the subsystem to (siSj)s and (Utj)t. In consequence we can distinguish \{ij)k > and 
\(ji)k > in this notation. In the partial wave decomposition there is a simple phase relation 
connecting both states 

\(ij)k >= n l+s+t urn >= r W) ik*)* > (A5) 

Of course all sets of basis states are complete, therefore we can expand the incoming state 
in several ways 

\ X > = ^|(31)2 ><(31)2|x> 

= ^|(12)3 ><(12)3|x> (A6) 

Due to the total antisymmetry of \x >, its matrix elements are equal 

< (31)2| X >=< (12)3| X >=< (23)1| X > (A7) 

In consequence we can expand the different terms in Eq. (|A2|) on the right hand side in 
different Jacobi coordinates. Let us begin with the two U terms in Eq.([A2|): 

M 3 = l - < (23)1|(12)3" >< (12)3"|T2(r 12 )|(12)3 w > 

< (12)3'"|(31)2"" > < (31)2' / ' / |T^(r 1 3)|(31)2 / > 
+ l - < (23)1|(31)2" > < (31)2"|T2(r 13 )|(31)2"' > 

< (31)2'"|(12)3"" > < (12)3""|T2(r 12 )|(12)3 / > (A8) 

In this equation we introduced several completeness relations and projected on two kinds of 
incoming states keeping in mind that one can add up the matrix elements because of the total 
antisymmetry of |% >. In the way we inserted the completeness relations in Eq . (|A8|) , the 
matrix elements of T 2 are evaluated in their natural coordinates. In this form the operator 
is diagonal in the quantum numbers and momenta of the outer particle and additionally it 
conserves the symmetry with respect to the interchange of the particles of the subsystem. 

The matrix element < (12)3|T^(r 12 )|(12)3' > depends on Jacobi coordinates and quan- 
tum numbers which single out the subsystem (12). By renumbering the particles one finds 

< (12)3|7; 2 (r 12 )|(12)3' >=< (31)2 |T^(r 13 ) | (31)2' > (A9) 

Note that the matrix element on the right hand side depends on coordinates which differ 
from the ones on the left hand side. Though the matrix elements are equal, one has to keep 
in mind that the actual meaning of the momenta and quantum numbers is different on both 
sides. 

In the same manner one can find an important relation between the cyclic and anti cyclic 
transformations 

< (ij)k\(jk)i' >=< iik)j\{kj)i' >= (_)(«)(-)<*>' < (ki)j\(jk)i' > (A10) 

We would like to emphasize that the transformation itself doesn't conserve the symmetry of 
the subsystem and therefore the completeness relation in " and "' states in Eq . ( |A8|) have to 
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include also the unphysical symmetric states in the two-body subsystems. With the help of 
Eq. ( |A10| ) and renumbering the particles in the second part of Eq. ( |A8| ) one finds 

M 3 = ~ < (23)1|(12)3" > < (12)3 /, |T w 2 (r 12 )|(12)3 /,/ > 

< (12)3'"|(31)2"" > < (31)2""|T 2 (r 13 )|(31)2' > 

+ I(_)(23)(_)(i2)"(_)(i2)'"(_)(3i)"" < (23)1|(12)3" > < (12)3"|i; 2 (r 12 )|(12)3 ,// > 

< (12)3" / |(31)2"" > < (31)2""|T 2 (r 13 )|(31)2' > (All) 

Because T 2 conserves the symmetry of the subsystem (— )( 12 )" = (— )( 12 )"' an d (— )( 31 )"" = 
(_)(3i) _ Therefore Eq. flAllj ) reduces to 

M 3 = \ (l + (-)( 23 )(-)( 31 )') < (23)1|(12)3" > < (12)3"|^(r 12 )|(12)3 ,// > 

< (12)3"'|(31)2"" > < (31)2""|T 2 (r 13 )|(31)2' > (A12) 

The incoming state is antisymmetric in the subsystem (31) hence (— )( 31 )' = —1. Therefore 
M 3 is zero for outgoing states which are symmetric in the subsystem (23) and equals twice 
the first part for the antisymmetric outgoing states. We restrict the outgoing states to 
antisymmetric ones. Then it is justified to write 

M 3 = < (23)1|(12)3" > < (12)3"|T 2 (r 12 )|(12)3'" > 

< (12)3"'|(31)2" ,/ > < (31)2""|T 2 (r 13 )|(31)2' > (A13) 

We would like emphasize again that we restrict the outgoing, incoming and "" states to 
physical, antisymmetric states in the subsystems. Due to the coordinate transformations, 
the " and '" sums are not restricted anymore to antisymmetric subsystem states and the 
completeness relations run over all symmetries. 

Lets turn to the A 2n parts now. In the same manner we define the matrix elements of 
the A 2n parts 

M 3 = < (23)1 1 (12)3" > < (12)3"|X 12 |(12)3 / " > 



< (12)3"'|2 (f 2 .r 3 + -f[-r 2 x r 3 j |(31)2"" > < (31)2""|X 13 |(31)2' > 
+ < (23)1 1 (31)2" > < (31)2"|X 13 |(31)2" / > 

< (31)2'"|2 (V 2 • r 3 - i H ■ r 2 x r{j |(12)3"" > < (12)3""|X 12 |(12)3 / > (A14) 

The isospin operators commute with the operators. We combined them with the inner 
transformation matrix. Because f 2 -f^, is symmetric with respect to the interchange of particle 
2 and 3 and f[ • f 2 X is antisymmetric, the steps leading to Eq. ([A10| ) can be repeated with 
the isospin matrix element. In addition there is a sign change for the antisymmetric part of 
the isospin operators. 

< (12)3"'|2 (rl ■ r 3 + % - t 1 ■ r 2 x r 3 ) |(31)2"" > 

\(31)"7 \(12)"" "'- f - - 1 



< (31)2"'|2 [f 2 ■ f 3 - - n ■ r 2 x r 3 ) |(12)3"" > (A15) 
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The symmetry of the Xi/s leads to equal phases on both sides of their matrix elements: 
( — )( 31 ) = (— )( 31 ) and (— )*- 12 ^ = (— ) ( - 12 - ) • The matrix element of the A 27r parts reduces to 

M 3 = (l + (_)(23)(_)(3in < (23)1|(12)3" > < (12)3"|X 12 |(12)3"' > 



< (12)3'"|2 (f 2 ■ r 3 + - f[ ■ r 2 x rlj |(31)2' w > < (31)2""|X 13 |(31)2 / > (A16) 

and its is again justified to restrict the outgoing states to antisymmetric subsystems and 
write 

M 3 = 2 < (23)1|(12)3" > < (12)3"|X 12 |(12)3" / > 



< (12)3"'|2 [r 2 ■ T3 + - f x ■ r 2 x r 3 j |(31)2"" > < (31)2""|X 13 |(31)2' > (A17) 

It is standard to work out the two-body partial wave matrix elements for Xy and T 2 in 
terms of Bessel transforms. Also the isospin and coordinate transformation matrix elements 



are standard and we refer to 48 
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TABLES 



A [m n ] 

CD Bonn+TM 4.856 
AV18+TM 5.215 
Nijm I+TM 5.120 
Nijm II+TM 5.072 
Nijm'93+TM 5.212 



TABLE I. The cut-off parameters A from Eq. (^) used in the given potential combinations. 
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FIGURES 




FIG. 1. The differential cross section in elastic Nd scattering at 3 (a), 65 (b), 135 (c) and 
190 MeV (d). Two bands are shown in each subfigure, the light shaded one contains NN force 
predictions, the darker shaded one the NN force predictions+TM 3NF. The solid curves are the 
AV18+URBANA IX predictions. The dashed curves are the CD Bonn+TM' predictions. Data 
at 3 MeV from |k| (pd), at 65 MeV from |6lJ (pd-crosses) and Q (nd-circles) , at 135 MeV 
from p2| (pd-crosses), [^3| (pd-circles), and at 190 MeV from Q (pd: crosses 181 MeV, circles 
216.5 MeV). In some cases error bars are not visible on the scale of the figure. 



17 




FIG. 2. The analyzing power A y {N) for elastic Nd scattering. Curves and the sequence of 
energies as in Fig. 0. Data at 3 MeV from |l| (pd -crosses) and f^4|] (nd-circles), at 65 MeV 
from |l| (pd-crosses) and || (nd-circles), at 135 MeV from |7| (pd-circles 150 MeV), @ 
(pd-crosses 146 MeV), [|| (pd-x's 155 MeV), and at 190 MeV from |27) (pd-crosses 190 MeV), 
|| (pd-circles 198 MeV), @ (pd-squares 197 MeV). 
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FIG. 3. The deuteron vector analyzing power iT%% for elastic Nd scattering. Curves and the 
sequence of energies as in Fig. |]. pd data at 3 MeV from [61], at 65 MeV from ]6^| , at 135 MeV 
from p2l (crosses), [E3] (circles), and at 190 MeV from 26]. 
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FIG. 4. The tensor analyzing power T20 for elastic Nd scattering. Curves and the sequence of 
energies as in Fig. |l]. pd data at 3 MeV from at 65 MeV from [|68| , and at 135 MeV from [22] 
(crosses), |23j (circles). 
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FIG. 5. The tensor analyzing power T21 for elastic Nd scattering. Curves and the sequence of 
energies as in Fig. |l]. pd data at 3 MeV from at 65 MeV from [|68| , and at 135 MeV from [22] 
(crosses), p3j (circles). 
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FIG. 6. The tensor analyzing power T22 for elastic Nd scattering. Curves and the sequence of 
energies as in Fig. |l]. pd data at 3 MeV from at 65 MeV from [|68| , and at 135 MeV from [22] 
(crosses), p3j (circles). 
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FIG. 8. The spin transfer coefficient K>j' y for elastic Nd scattering. Curves and the sequence of 
energies as in Fig. |]. pd data at 135 MeV from [f2^]. 
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FIG. 9. The spin transfer coefficient K^.' x for elastic Nd scattering. Curves and the sequence of 
energies as in Fig. |]. pd data at 135 MeV from 
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FIG. 10. The spin transfer coefficient K^ z for elastic Nd scattering. Curves and the sequence 
of energies as in Fig. |l[ pd data at 135 MeV from [22]. 
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FIG. 11. The spin transfer coefficient K^ z for elastic Nd scattering. Curves and the sequence 
of energies as in Fig. |l[ 



27 




30 60 90 120 150 180 30 60 90 120 150 180 



9 [deg] 9 [deg] 



FIG. 12. The spin correlation coefficient C XVjX for elastic Nd scattering. Curves and the se- 
quence of energies as in Fig. |l[ 
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FIG. 13. The spin correlation coefficient C yy for elastic Nd scattering. Curves and the sequence 
of energies as in Fig. |l[ pd data at 190 MeV from [26]. 
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FIG. 14. The spin correlation coefficient C zz for elastic Nd scattering. Curves and the sequence 
of energies as in Fig. |l[ 
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